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Factorial designs have broad applications in agricultural, engi¬ 
neering and scientific studies. In constructing and studying proper¬ 
ties of factorial designs, traditional design theory treats all factors 
as nominal. However, this is not appropriate for experiments that in¬ 
volve quantitative factors. For designs with quantitative factors, level 
permutation of one or more factors in a design matrix could result 
in different geometric structures, and, thus, different design proper¬ 
ties. In this paper indicator functions are introduced to represent 
factorial designs. A polynomial form of indicator functions is used to 
characterize the geometric structure of those designs. Geometric iso¬ 
morphism is defined for classifying designs with quantitative factors. 

Based on indicator functions, a new aberration criteria is proposed 
and some minimum aberration designs are presented. 

1. Introduction. Factorial designs are commonly used in most industrial 
and scientific studies. In such a study, a number of fixed levels (settings) are 
selected for each factor (variable), and then some level combinations are 
chosen to be the runs in an experiment. A factor can be either nominal or 
quantitative. For nominal factors, there is no ordering among levels. The 
interest of analysis of an experiment with nominal factors is to understand 
if there exist differences in treatment means and if they exist, which treat¬ 
ment means differ. Analysis such as ANOVA or various multiple comparison 
testing procedures is often used for treatment comparison. In many studies, 
especially in response surface exploration, factors are often quantitative and 
there exists an order among levels. For an experiment with quantitative fac¬ 
tors, the objective is usually achieved through fitting a (polynomial) model 
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Table 1 

Combinatorially isomorphic 
designs with different geometric 
structures 
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that can “well” describe the relationship between the response and the fac¬ 
tors. The distinction in the analysis objective and strategy for these two 
types of experiments requires different selection criteria and classification 
methods. 

For designs with nominal factors, the design properties should be invariant 
to level permutation within one or more of its factors. However, for quantita¬ 
tive factors, Cheng and Wu (2001) observed that level permutation of 3^“^ 
designs could result in changes in model efficiency when a polynomial model 
is htted, which is referred to as “model nonisomorphism.” Independently, Ye 
(1999) also observed that level permutation could alter the aliasing structure 
of designs when linear-quadratic decomposition [see Wu and Hamada (2000), 



Fig. 1. Combinatorially isomorphic but geometrically nonisomorphic designs. 
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Section 5.6] is used. As shown in the following example, such “model noniso¬ 
morphism” is, indeed, the result of different geometric structures induced by 
permuting levels of factors. Consider the two 3^“^ designs in Table 1. In the 
table each design is written as a design matrix in which each column repre¬ 
sents a factor and each row represents an experimental run. These designs 
are combinatorially isomorphic since one is obtained by applying the permu¬ 
tation {0,1,2} —> {0,2,1} on the third column of the other. However, if we 
treat these levels as quantitative, their geometric structures are apparently 
different as shown in Figure 1. The difference in geometric structure also 
reflects on the model efficiency. For example, when a model that contains 
all linear main effects and three linear-by-linear interactions is considered, 
the design on the left-hand side has higher D-efficiency than the one on the 
right-hand side. 

The conventional mathematical tools used for factorial designs, such as 
group theory and coding theory, treat all factors as nominal. Therefore, they 
do not differentiate geometric structures resulting from level permutations 
and fail to study the design properties associated with its geometric struc¬ 
ture. 

A new approach for characterizing designs with quantitative factors is de¬ 
veloped in this paper. When all k factors in a factorial design are quantita¬ 
tive, it can be viewed as a collection of points in M^. This collection of points 
is represented by an indicator function, which will be defined in Section 2. 
The indicator function can be written in a polynomial form which reveals 
the design’s properties and characterizes its geometric structure. Thus, clas- 
sihcation and design criteria are developed based on the indicator functions. 
This approach is motivated by Pistone and Wynn (1996), which first used 
polynomial systems to describe designs and studied their properties using 
algebraic geometry methods. In this paper properties of designs with quan¬ 
titative factors are studied. Section 2 introduces the indicator function as 
a mathematical tool for examining the geometric structures of designs. In 
Section 3 geometric isomorphism is dehned for the classification of factorial 
designs. Section 4 proposes a new aberration criterion for factorial designs 
with quantitative factors. Some remarks are given in Section 5. 

In the remainder of this section we will introduce some notation and 
terminology. Let T> be the OA(A^, siS 2 ■ • • Sk), which is a full factorial design 
with k factors and N design points, where N = siS 2 ■ ■ - Sk- Unless specified, 
the levels of ith factor are set at Gj = {0,1,..., Sj — 1} C M for each factor, 
which are evenly spaced. Therefore, 2? is a set of N points in A fe-factor 
factorial design A is said to be in a design space V if its design points are 
all in T), that is, Vx G A, x G P. A design point in T) may appear more than 
once in A. Throughout this paper, X]xeA/(^) sums the function / over all 
design points in A] that is, if x appears multiple times, /(x) is summed over 
multiple times. 
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For each factor Xi, define a set of orthogonal contrasts Cq{x),CI{x), ..., Cl._ 
such that 


( 1 . 1 ) 


^ Clix)Cl{x) 


0, if V, 

Si, if u = v. 


Let T = Gi X • • • X Gfc. An orthonormal contrast basis (OCB) on T> is defined 
as 

k 

(1.2) Ct(x) = nQ^(x,) 

i=l 


for t = {ti,t 2 ,... ,tk) gT and x = {xi,X 2 , ■ ■ ■, Xk) G V. It is obvious that 


(1.3) 


Ct(x)Cu(x) 


0, if t / u, 

N, if t = u, 


where t, u are elements in the set T. In statistical analysis, C'q(x) = 1 is often 
adopted to represent a constant term. Therefore, we call {C't(x)} with Cq = 1 
for all i a statistical orthonormal contrast basis (SOCB). When Cj{x) is a 
polynomial of degree j for j = 0,1,..., Sj — 1 and i = 1,2,..., fc, the SOCB 
is called an orthogonal polynomial basis (OPB) [Draper and Smith (1998), 
Chapter 22]. Note that an OPB is an SOCB, and an SOCB is an OCB. 

We define two norms on T. Let ||t||o be the number of nonzero elements 
in t and let 

k 

i=l 

For a contrast Ct in an SOCB, ||t||o is the number of factors it involves. If 
the SOCB is also an OPB, ||t||i gives its polynomial degree. Two contrasts 
Ct and Cu in an SOCB are said to be disjoint if they have no common 
factors, that is, maxi<j<fcmin(ti,Uj) =0. 


2. Indicator functions. Indicator functions are presented in Fontana, Pi- 
stone and Rogantin (2000) for studying two-level fractional factorial designs 
(without replicates). Ye (2003) generalizes to accommodate replicates. In 
this section the definition is extended further to general factorial designs. 

Definition 2.1. Let A be a design in the design space D. The indicator 
function of A is a function defined on D, such that for x gT>, the 

value of is the number of appearances of point x in design A. 

The following proposition follows immediately from the definition. 
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Proposition 2.1. LetAi,...,Am be factorial designs of the same de¬ 
sign space V and Fyi. (x), i = 1, ... ,m, be their corresponding indicator func¬ 
tions. Let B be the combined design (design points are repeatable in B) of 
All,..., Am- Then the indicator function of B is 

m 

^b(x) =^PU,(x). 
i=l 

Since a design is uniquely represented by its indicator function, the indi¬ 
cator function carries all properties of this design. Some of these properties 
are revealed when indicator functions are expanded with respect to an OCB. 

Theorem 2.1. Let A be a faetorial design with n runs. Let V be the de¬ 
sign spaee of A, and {C't(x),t G T} be an OCB defined on V. The indieator 
funetion of A ean be represented as a linear eombination of Cts as follows: 

(2.1) F^(x) = ^6tCt(x), 

t&T 

for all X G 2?. The eoeffieients {5t,t G T} are uniquely determined as 

(2.2) ^t = ^E^t(x). 

And, in particular, for an SOCB, bo = n/N, where 0 = (0,0,..., 0). 

Proof. The indicator function is defined on B and Fj[(T>) can 

be viewed as a vector in M'^. Since the {Ct{T)),t G T} forms a basis of M^, 
Fj[{'D) can be represented as a linear combination of {Ct{T))'\. Equivalently, 
(2.1) is true. For the coefficients bts, 

^ Ct(x) = F^(x)C-t(x) = E E &sCs(x)Ct(x) 

xGa4 x£X) xG(D s£T 

= E^sECs(x)C't(x) = iV6t. 
ssT xev 

The proof is complete. □ 

Note that the theorem does not depend on level settings and choice of 
Cj{x), as long as (1.1) is satisfied and {Ct} is defined as in (1.2). In the 
functional space generated by linear combinations of {Ct}, the indicator 
function of a design has a unique representation, that is, there is a one-to- 
one relation between a factorial design and its bt values. This is an extension 
of a similar result on two-level designs presented in Ye (2003). When {Ct} is 
an OPB, an indicator function can be uniquely represented as a polynomial 
of degree no more than nf=i('S* “ 1). 
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A projected design has the same number of runs as the original design but 
is in a reduced design space with only a subset of the original factors. Given 
a design’s polynomial representation in the form of (2.1), the polynomial 
representations of its projected designs are easily available, as shown in the 
following corollary. 

Corollary 2.1. Let A be a factorial design in design space T) and 
Fyx{x.) = ^tC't(x) be its indicator function. Without loss of generality, 
let B he its projection to factors Xi,... ,Xi. If {Ct} is an SOCB, the indi¬ 
cator function of B is then 

(2.3) FB{xi,...,xi)=N2'^btCt, 


^2= n 'Fi = {t\ti+i = ■ ■ ■ = tk = 0}. 

i=l+l 

Proof. From (2.2), bt = ^t(x). coefficient of C't(x) in 

Fg(x) is then 1/A^i ^t(x), where A^i = 0^=1 Equation (2.3) follows. 

□ 


The coefficients bt also relate to the orthogonality of a design. This can 
be shown in the following corollary which follows immediately from (1.2) 
and (2.2). 

Corollary 2.2. Let {C't(x),t sT} be an SOCB. For disjoint Cu and 

Cv, 

Wv = ^ C’u(x)Cv(x). 

' xeA 

Furthermore, the correlation of Cu and in A is 6u+v/^o- 

Let C'u(x) and C'v(x) be two disjoint contrasts. From Corollary 2.2, the 
two contrasts are zero correlated on design A if and only if 6u+v = 0. As a 
special case, bt = 0 implies that the contrast C't(x) has zero correlation with 
the constant term on design A. In general, a smaller bt implies a lesser degree 
of aliasing between effects and, therefore, bt can be used as a measurement 
of aliasing between effects. Various statistical properties of designs can be 
studied through the 6t’s. More results will be shown in Sections 3 and 4. 
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Example 2.1. Consider the case of three-level factorial designs with 
k factors. The design space V is a collection of 3^ points: {(di,..., 
dj = 0, 1,2, z = 1,... , k}. From Definition 2.1, any fc-factor three-level facto¬ 
rial design can be represented by an indicator function defined on D. The 
orthonormal polynomials for a three-level factor are 

Co(a;) = 1, Ci{x) = \J\{x - 1) and C2(x) = \/2(|(x - 1)^ - 1). 

Note that (Ci(0), 6-1(1), Ci(2)) = and (62(0), 6-2(1), C2(2)) 

(l/\/2, —\/2, l/\/2) are proportional to the linear and quadratic contrasts, 
respectively, as dehned in Wu and Hamada [(2000), Section 5.6]. Thus, 
{6t(x),t G T}, where T is the vector space {0,1,2}^, is an OPB for the 
functional space of T). By Theorem 2.1, an indicator function can be written 
as a linear combination of 6t(x)’s with coefficients 

xeA 

and, in particular, bo = nj"^. The coefficients bt contain information about 
aliasing between effects. For the design on the right-hand side of Table 1, 
which is also shown on the right-hand side of Figure 1, its indicator function 
is 

F(x) = i6ooo(x) - ^6111 (x) - ^6 ii2(x) + ^6121 (x) + ^6-2ii(x) 

- ^Ci22(x) - ^C212(x) + ^6221 (x) + ^6222 (x), 

where 


62122^3 (x) (^1) ^ i 2 (^2)623 (X3). 

For the design on the left-hand side, its indicator function is 

F(x) = i6ooo(x) + ^6 ii2(x) + ^Ci2l(x) + ^C21l(x) - ^6222 (x). 

(2.5) 

In (2.5), 6111 = 0 implies that the linear-by-linear interactions are orthogo¬ 
nal to linear main effects in the design. For the other design, they are not 
orthogonal since bm = -^ in (2.4). This is consistent with the higher D- 
efficiency of the former design when a model with all linear main effects and 
linear-by-linear interactions is considered. 

3. Geometric isomorphism. When factors are all nominal in a factorial 
design, new design matrices obtained through level permutations in one or 
more factors are considered to be isomorphic to the original design. This is 
referred to as the combinatorial isomorphism. As shown in Section 1, when 
factors are quantitative, level permutations generate designs with different 
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geometric structures and, thus, different design properties. Cheng and Wu 
(2001) observed differences in H-efficiency of these designs and proposed 
model isomorphism for classihcation. However, such classification depends 
on a priori specihed models. Designs that have the same efficiencies with 
respect to a given model might have different efficiencies with respect to 
another model. A classification with respect to a certain model can be no 
longer useful when a different model is considered. For consistency, a better 
classihcation method should be based on the geometric structures, which 
are fundamental to design properties and do not depend on the choice of 
the models. 

From a geometric viewpoint, a geometric object remains the same struc¬ 
ture when rotating and/or rehecting with respect to a super-plane. In the 
context of a factorial design, only rotations and rehections, after which the 
resulting designs are still in the design space V, should be considered. Ro¬ 
tating then corresponds to variable exchange and rehecting corresponds to 
reversing order of the levels. Therefore, we dehne geometric isomorphism of 
two designs as follows. 

Definition 3.1. Let A and B be two factorial designs from the same 
design space V. Designs A and B are said to be geometrically isomorphic if 
one can be obtained from the other by variable exchange and/or reversing 
the level order of one or more factors. 

One can differentiate geometrically nonisomorphic designs by comparing 
their indicator functions. Let ... ,Xk) be the indicator function of 

design A and Ai be the design obtained by reversing the level order of factor 
Xi in A. The indicator function of Ai is Fj\^{2d — xi,X 2 , ■ ■ ■ ,Xk), where d is 
the center of the levels. Let A 2 be the design obtained by exchanging factor 
Xi with factor X 2 in A. Then its indicator function is Fj[(x 2 ,xi,... ,Xk)- 
If the indicator functions of two designs are the same after a series of such 
operations, then they are geometrically isomorphic. 

Geometric isomorphism of two designs can be more easily examined when 
indicator functions are expanded to polynomial form with respect to an 
OPB. Theorem 3.1 implies that if two designs are geometrically isomorphic, 
the absolute values of their coefficients bt must show the same frequency 
patterns. 

Theorem 3.1. Let A and B be two factorial designs of the design 
space V, and {C't(x)} he an OPB defined on V. Let F_a{-x) =X)^tC't(x) 
and Fb(x) =X)^tC't(x) be the indicator functions of A and B, respectively. 
Designs A and B are geometrically isomorphic if and only if there exist a 
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permutation {iii 2 ■ ■ - ik) and a vector {jij 2 ■ ■ -jk), where ji’s are either 0 or 
1, such that 



(3.1) 


for all t = {tit 2 .. .tk) G T. 

Proof. Using the three-term recursive equation given in Kennedy and 
Gentle [(1980), pages 343 and 344] for constructing orthogonal polynomials, 
it is easy to show that the orthogonal polynomial contrasts Cj{x) of a factor 
satisfy the following condition: 



if j odd, 
if j even. 


(3.2) 


If A and B are geometrically isomorphic, then by definition A must be ob¬ 
tained from B by variable exchange and/or reversal of levels. Let the variable 
exchange be xi —> Xj,, and let ji = 1 if the levels of factor xi are reversed, 
ji = 0 if not. Hence, (3.1) is truly based on (2.2) and (3.2). Conversely, if 

(3.1) is true, B can be obtained from A by the variable exchange xi —> Xj; 
and the level reverses on the factors with ji = 1. Therefore, A and B are 
geometrically isomorphic. □ 

Note that from the proof Theorem 3.1, it holds for any basis such that 

(3.2) is satished. From the theorem, one can immediately show that (2.4) 
and (2.5) represent two geometrically nonisomorphic designs as their coeffi¬ 
cients show different frequency patterns. In general, with a proper choice of 
{Ct(x)}, two designs are geometrically isomorphic if and only if their indi¬ 
cator functions have the same coefficients bt after a certain type of permu¬ 
tation and sign reversal. Otherwise, if two designs have different geometric 
structures, their coefficients must show different frequency patterns. 

Example 3.1. The Lig array (shown in Table 2) is one of the most 
popular designs among industrial experimenters. It has one two-level factor 
and seven three-level factors. For the moment, we only consider the three- 
level factors. Wang and Wu (1995) studied the projected design of Lis and 
reported three combinatorially nonisomorphic cases for 3-factor projections 
(denoted as 18-3.1, 18-3.2 and 18-3.3) and four combinatorially nonisomor¬ 
phic cases for 4-factor projections (denoted as 18-4.1, 18-4.2, 18-4.3, 18-4.4). 
As shown earlier in this paper, level permutation in a design may create 
designs with different geometric structures. There are a total of six permu¬ 
tations among three levels, that is, 
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The six permutations can be divided into three pairs as shown in Table 
3. Within each pair, one permutation is the reverse of the other, hence, 
only one is needed in generating geometrically nonisomorphic designs. For 
each combinatorially nonisomorphic case, permutations are applied to each 
column to search for all geometrically nonisomorphic cases. For three-factor 
projections, there are two, four and two geometrically nonisomorphic cases 
within 18-3.1, 18-3.2 and 18-3.3, respectively. Cheng and Wu (2001) reported 
the same number of model nonisomorphic cases for 18-3.1 and 18-3.2 but did 
not report model nonisomorphic cases for 18-3.3. For four-factor projection, 
there are four, ten, three and four geometrically nonisomorphic cases in 18.4- 
1, 18.4-2, 18-4.3 and 18-4.4, respectively. Cheng and Wu (2001) only reported 
four model nonisomorphic cases for 18-4.2 and none for the other three. A 
complete list of these geometrically nonisomorphic projected designs is given 
in the Appendix. 


4. Aberration criterion. A popular criteria for factorial designs is min¬ 
imum aberration. The original definition of minimum aberration based on 
group theory applies to regular p'^~^ fractional factorial designs [Fries and 
Hunter (1980)]. Recently, Xu and Wu (2001) proposed an aberration crite¬ 
rion based on coding theory for general factorial designs. It reduces to the 


Table 2 

Lis orthogonal array 


0 

1 

2 

3 

4 

5 

6 

7 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

1 

1 

1 

1 

1 

1 

0 

0 

2 

2 

2 

2 

2 

2 

0 

1 

0 

0 

1 

1 

2 

2 

0 

1 

1 

1 

2 

2 

0 

0 

0 

1 

2 

2 

0 

0 

1 

1 

0 

2 

0 

1 

0 

2 

1 

2 

0 

2 

1 

2 

1 

0 

2 

0 

0 

2 

2 

0 

2 

1 

0 

1 

1 

0 

0 

2 

2 

1 

1 

0 

1 

0 

1 

0 

0 

2 

2 

1 

1 

0 

2 

1 

1 

0 

0 

2 

1 

1 

0 

1 

2 

0 

2 

1 

1 

1 

1 

2 

0 

1 

0 

2 

1 

1 

2 

0 

1 

2 

1 

0 

1 

2 

0 

2 

1 

2 

0 

1 

1 

2 

1 

0 

2 

0 

1 

2 

1 

2 

2 

1 

0 

1 

2 

0 
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Table 3 
Six 

permutations of 
three levels 


I 

u 


0 

I 

2 

I 

2 

0 

2 

0 

I 

rv 

rv 

rv 

2 

I 

0 

I 

0 

2 

0 

2 

I 


traditional aberration criterion for regular p^~'^ designs, and the G 2 aberra¬ 
tion criteria [Tang and Deng (1999)] for general two-level factorial designs. A 
statistical justification is given by Xu and Wu (2001) to relate the criterion 
with ANOVA. From this relation, it can be easily seen that their aberration 
criterion can be redefined using the indicator functions as follows. 

Definition 4.1. Let ^ be an n x A: factorial design of design space 
V. Let Fyx{x) = ^tC't(x) be its indicator function, where {Ct} is an 
SOCB. The generalized wordlength pattern (ai(.A),... ,ak{A)) of design A 
is dehned as 

(4.1) «.o)=T,(^)'- 

II t II0—^ 

The generalized minimum aberration criterion is to sequentially minimize 
ai{A) for i = 1,2,..., A:. The resolution of A equals the smallest r such that 
ar > 0 . 

From Corollary 2.2, (6t/^o)^ is a measurement that reflects the sever¬ 
ity of aliasing between the effect Ct and the general mean. Therefore, in 

(4.1) , Oi measures the overall aliasing between all f-factor effects and the 
general mean. A smaller ai indicates a lesser degree of aliasing between the 
z-factor effects and the overall mean. Therefore, the a^’s are to be minimized 
sequentially. Note that the definition implicitly assumes all z-factor effects 
are equally important, which is only suitable for nominal factors (see later 
discussion on the hierarchical ordering principle). The above definition is a 
natural generalization of the definition of the aberration criterion for two- 
level factorial designs given in Ye (2003). Xu and Wu (2001) showed that 
A is an orthogonal array of strength t if and only if ai{A) = 0 for 1 < z < t. 
From the definition, the if and only if condition can be stated in the language 
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of indicator function as follows: 6t = 0 for all t’s such that 1 < ||t||o < t- For 
example, the coefficients 6t in indicator functions (2.4) and (2.5) are zero 
for all t such that 1 < ||t||o < 2. Therefore, they are orthogonal arrays of 
strength two. This aberration criterion is invariant to level permutation as 
well as the choice of contrasts. While these features are desirable when all 
factors are nominal, it is not quite desirable when factors are quantitative. 
An immediate problem is that this aberration criterion completely fails to 
distinguish and rank combinatorially isomorphic but geometrically noniso¬ 
morphic designs, for example, the two designs in Table 1. 

An important assumption behind aberration criteria is the hierarchical 
ordering principle [Wu and Hamada (2000), Section 3.5]: (i) low-order effects 
are more likely to be important than high-order effects, and (ii) effects of the 
same order are equally likely to be important. The principle can be applied to 
nominal and quantitative factors. However, the effect orders for the two types 
of factors should be different. For nominal factors, the objective of analysis is 
treatment comparison. Therefore, all i-factor effects are regarded as equally 
important and z-factor effects are more important than j-factor effects for 
i <j. The effect order is decided by the number of factors that are related 
to the corresponding contrast Ct, that is, the value of ||t||o. Therefore, in 

(4.1) , the overall aliasing is measured by taking the sum over those t’s with 
the same ||t||o value. For experiments with quantitative factors, polynomial 
models are often utilized to approximate the response. In this case, effects 
of higher polynomial degree are regarded as less important than effects of 
lower polynomial degree. Therefore, the order of effect importance should 
be arranged according to polynomial degrees. Recall that in an OPB, Cj{x) 
is a polynomial of degree j and C't(x) is a polynomial of degree ||t||i. In this 
case, the order of effect importance can be defined according to the values 
of ||t||i. For example, for quantitative three-level factors, the order of effect 
importance is as follows: 

I » q == ll » Iq == ql == III 

(4.2) 

>> qq == llq == Iql == qll == llll » • • •, 

where >> is read as “more important than” and == as “as important 
as,” and I and q indicate linear and quadratic main effects, respectively, 
ll linear-by-linear interaction, etc. Such ordering is consistent with the re¬ 
sponse surface methodology in which, based on a rationale from Taylor’s 
series expansion, effects of the same degree are sequentially added to the 
model starting from the lowest degree. 

Consider the two designs in Table I. When all factors are quantitative and 
OPB is used, the contrasts in equations (2.4) and (2.5) follow the linear- 
quadratic decomposition and have clear interpretation in terms of htting 
polynomial models. As mentioned previously, in (2.5), = 0 implies that 
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the III interactions have zero correlation with the constant term, and the ll 
interactions between any two factors have zero correlation with the linear 
main effect of the third factor. However, the design of (2.4) does not have 
this nice property since 6111 = hence, the III interaction is (partially) 
aliased with the constant term, similar with the linear main effects and 
ll interactions. In Xu and Wu (2001), hm is considered as important as 
^ 112 , ■ • ■, ^222 and are summed up together in as. This is not quite appropriate 
if a polynomial model is to be fitted and III interactions are regarded as more 
important than all other three-way interactions as in (4.2). Therefore, when 
all factors are quantitative with polynomial models as our point of interest, 
the following aberration and resolution criteria are proposed. 

Definition 4.2. Let ^ be an n x fc factorial design with quantitative 
factors, and let {Ct} be an OPB. Let T 4 (x) = ^tC't(x) be the indicator 
function of A. The generalized wordlength pattern {f3i{A),..., Pk{A)) is 
defined as 

(4.3) f3,{A) = 

||t||i =2 

The generalized minimum aberration criterion is to sequentially minimize 
f3i for f = 1 , 2 ,..., LC, where 

k 

K = '£{si-l) 
i=l 

is the highest possible degree in the decomposition. The resolution of A is 
defined to be the smallest r such that (Sr > 0 . 

In the rest of this paper we refer to the word length pattern given in 
Definition 4.1 as the a wordlength pattern, and the wordlength pattern 
in the above definition as the (5 wordlength pattern. It is probably more 
appropriate to call the new aberration and resolution criteria polynomial 
degree” aberration and resolution. For three-level designs, (4.3) counts the 
overall aliasing between the general mean and effects that are of the same 
importance in (4.2). Based on the above definition, the wordlength patterns 
of the two designs in (2.4) and (2.5) are (0,0, |, |, |) and (0,0,0, |,0, ^), 

respectively. The latter one has less aberration and higher resolution and is 
favored. 

Two contrasts, C't(x) and C'u(x) G T, are treated as equally important if 
and only if ||t||i = ||u||i. Therefore, by Theorem 3.1, we can easily show that 
geometrically isomorphic designs have identical j3 wordlength patterns. 

Corollary 4.1. Let A and B be two geometrically isomorphic designs 
in design space D. Then their f5 wordlength patterns are identical. 
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Proof. Let the indicator functions of two designs be = '^htCt{x) 
and Fq = J2 By Theorem 3.1, there must be a permutation (iiZ 2 ■ ■ - ik) 

such that for all t = (tit 2 ■ ■ - tk) ^ T. For each t, 

denote its corresponding permutation as t'. Since ||t||i = ||t'||i, (5j{A) and 
PjiB) sum over the same values, and hence, are identical. □ 

By Corollary 4.1, two designs with different (5 wordlength patterns must 
be geometrically nonisomorphic. 

The only distinction between Definitions 4.1 and 4.2 are the norms of t 
used in (4.1) and (4.3), which reflect the difference in ordering effects for 
nominal and quantitative factors. For a given design, the sum of its at’s is 
the same as the sum of its /3j’s. The following theorem shows that this sum is 
a constant for designs that have the same run sizes and replication patterns. 

Theorem 4.1. Let A be an n x k factorial design in the design space 
V. Let F_a{x) = X^teT ^tC't(x) he the indicator function of A. Then 



where 


n 2 ='^F^{x.) and N = si...Sk- 

xGX> 

For designs with no replicates, n 2 = n. 

Proof. 

^Fj(x) = ^('^6tCt) = E E &ti6t.Ct,(x)Ct,(x) 

xGD xgX> VtgT / xgDti,t 2 ST 

= E bt,bt2j2Ct,{^)Ct,{x) = Nj2bl 

ti,t 2 eT xex> teT 

From Theorem 2.1, bo = n/N. Hence, (4.4) is obtained. For designs with no 
replicates, ^^(x) = F 4 (x), hence n 2 = n. □ 

A special case of the above theorem is a well-known result for regular 
fractional factorial designs, in which the sum of their wordlength pat¬ 

tern vector equals p™ — 1. The theorem shows that it holds for both a and 
P wordlength patterns. The theorem also shows that the sum of wordlength 
pattern vectors is larger for designs with higher degrees of replication as n 2 
is larger in (4.4). Therefore, they tend to have higher aberration than those 
with less replicates. 
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Table 4 

Minimum aberration projections of Lis with only three-level factors 


ff of factors 

Columns 

{fia, fin) 

Resolution 

3 

W 2 5 

(0, 0.125, 0.75) 

IV 

4 

lu^ 2u^ 5 

(0, 1.875, 0) 

IV 

5 

lu^ 2u^ 3u^ 4 5 

(0, 6.0625, 0) 

IV 

6 

2u 3tt 4u 5 6 7 

(0.75, 6.9375, 6.75) 

III 

7 

In 2u 3u 4u 5 6 7 

(1.5, 14.625, 12) 

III 


Note that for two-level designs, both a and /3 wordlength patterns reduce 
to the same generalized wordlength pattern by Tang and Deng (1999) and 
Ye (2003). However, if the factors have more than two levels, these two 
wordlength patterns often give different minimum aberration designs. One 
should choose from them based on the nature of the factors, nominal or 
quantitative. 

Example 4.1. When less than seven three-level factors are considered 
in an experiment, it would be of interest to know which columns in the Lig 
array are the best to be assigned to those factors. To find the minimum 
aberration projections of the Lig array, an exhaustive search over all pos¬ 
sible projections was performed, and three-level permutations were applied 
to each column. Tables 4 and 5 list the (5 minimum aberration projected 
designs with and without the two-level factors, respectively. In the tables, 
u denotes the permutation {0,1,2} —> {1,2,0}; v? denotes the permutation 
(0,1,2} —> {2,0,1}. For example, the best projection with three 3-level fac¬ 
tors is columns Iv?, 2 and 5, where lu^ means permutation applies to 
column 1 in Table 2. It should be mentioned that, with exception of one 
trivial case in which the full factorial design is the only nonisomorphic de¬ 
sign, none of the designs in Tables 4 and 5 is combinatorially isomorphic to 
the minimum aberration designs given by Xu and Wu (2001). 


Table 5 

Minimum aberration projections of Lis with the two-level factors 


ff of factors 

Columns 

{fis,f3A, /3 b) 

Resolution 

3 

0 1 2 

(0, 0, 0) 


4 

0 1m^ 2 5 

(0, 0.5, 1) 

IV 

5 

0 1 3m 4 7 

(0, 3.75, 0) 

IV 

6 

0 1m 2 Sm^ 4m^ 7 

(0, 10.0625, 0) 

IV 

7 

0 1m 2 3u^ 4m^ 5 7 

(1.25, 14.21874, 7.40625) 

III 

8 

0 1 2m 3m 4 5m 6 7 

(2.5, 22.5, 17.3125) 

III 
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5. Concluding remarks and some discussion. This paper proposes a ge¬ 
ometric approach in studying factorial designs with quantitative factors. 
When factors are quantitative, the traditional mathematical treatment, which 
is appropriate for designs with nominal factors, no longer applies. The key in 
our approach is indicator functions and their polynomial forms as expanded 
to OCBs. They are used to distinguish designs’ geometric structures, which 
carry the designs’ properties. This approach is still appropriate even for 
ordinal categorical factors. 

In Corollary 2.2, the connection between bt and the aliasing of contrasts 
with no common factors is given. For contrasts that are not disjoint, the 
calculation of their aliasing (correlation) is more complex, but still depends 
on the 6t’s. In the following, a general formula is offered, which covers the 
situations of disjoint and nondisjoint contrasts. Let {Ct(x),t G T} be an 
SOCB. For each factor Xi, any product of its two contrasts C* (x)C'*(x) can 
be expressed as a linear combination of Cq{x), Cl(x), ..., and Cl._i{x) on 
the space {0,..., Sj — 1}. Let 

Cuix)Cl{x) = ^ for x = 0,1,...,Si - 1. 

w=0 

The correlation of two contrasts C'u(x) and C'v(x) can be written as a linear 
combination of 6t’s by the following formula: 


1 Cu(x)C-v(x) 
xeA 


(5.1) 




N 


Xi 


xeA *=i 


k Si —1 


ERE 


xSAj=l Wi=0 


si-1 S2-1 / k 

E E E n 

Wl=0W2=0 \2 = 1 





E <^w(x) 

xsA 


Note that for disjoint contrasts Cu and Cv, where u = (ui,..., Uk) and v = 
{vi,...,Vk), 


(5.2) 


Ui Vi, 

\ 0, otherwise. 


for i = 1,... ,k, and Corollary 2.2 can be also derived from equations (5.1) 
and (5.2). To demonstrate the calculation of aliasing between nondisjoint 
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contrasts using (5.1), let us consider the two designs in Example 2.1. Notice 
that for X € {0,1,2}, 

Ci{x)Ci{x) = j^C2{x) + l, 

Ci{x)C 2 {x) = ^Ci{x) and 
C2{x)C2ix) = ^C2{x) + l. 

Therefore, by (5.1) the correlation between two nondisjoint contrasts, say 
Ciio(x) andC'ioi(x), equals (;^ 62 ii+ 6 oii)/ 6 ooo- Its value is j for (2.4) and ^ 
for (2.5). In general, the aliasing among nondisjoint contrasts has also been 
captured in the wordlength patterns, which can be viewed as a summary 
measure of aliasing and are easy to compute. In theory, one can derive a 
criterion that explicitly calculates aliasing among all pairs of contrasts by 
laborious computation and this deserves some further investigation. 

The (5 wordlength pattern can be generalized when effect order is defined 
in other ways. For example, the effect order in Definition 4.1 is based on 
the number of factors that correspond to an effect, whereas the effect order 
in Definition 4.2 is based on the degree of the polynomial that represents 
an effect. The two characteristics can be combined to rank effect orders as 
follows: (a) first use the degree of polynomials to rank effects and then for 
those effects with the same order, use the number of factors to further rank 
their order; or (b) first use the number of factors to rank effects and then for 
those effects with the same order, use the degree of polynomials to further 
rank their order. For three-level designs, (a) generates the following order: 

I »q»ll »lq==ql 
>> III » qq » llq == Iql == qll » llll » • • •, 
and for (b), the effect order follows: 

I »q»ll »lq==ql 

>> qq >> III » llq == Iql == qll » llll » ■ ■ ■. 

In either case, the wordlength patterns can be defined by taking the sum 
of {bt/bo)'^ over coefficients of the Ct’s that are considered equally impor¬ 
tant. The corresponding aberration criteria sequentially minimize these sums 
starting from the most important effects. In general, this methodology is 
very flexible and can be applied on any reasonable effect orders. Note that 
Corollary 4.1 and Theorem 4.1 still hold under these wordlength patterns. 
In addition, although j3 wordlength patterns are defined on the coefficients 
with respect to an OPB, they can be defined with respect to an SOCB as 
long as an appropriate effect ordering exists. 

We chose to present this work in a self-contained fashion rather than 
with full algebraic geometry language, so that the ideas are more accessible 
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to the statistical community. Nonetheless, this work is another example of 
how algebraic geometry can be applied to statistics, and we will continue 
to explore the connections between the two helds. For other applications 
of algebraic geometry methods in statistics, see Pistone, Riccomagno and 
Wynn (2000). 


APPENDIX 

The geometrically nonisomorphic projected designs of Lig are listed in 
Table 6. All but two pairs of these designs have distinct wordlength patterns 
as dehned in (4.3). 


Table 6 

Combinatorially and geometrically nonisomorphic projected designs of 

Lis 


Comb, nonisomorphic 

Geom. nonisomorphic 

WLP (/33,/34,/36) 

18-3.1 

12 3 

(0.09375,0.09375,0.2813) 


lu^ 2 3 

(0,0.375,0) 

18-3.2 

12 5 

(0.09375,0.594,0.281) 


lu^ 2 5 

(0,0.125,0.75) 


lu 2u^ 5 

(0.375,0.125,0.375) 


lu^ 2u^ 5 

(0,0.5,0) 

18-3.3 

13 4 

(0.375,0.375,1.125) 


lu 3 4 

(0,1.5,0) 

18-4.1 

2 3 4 5 

(0.375,0.515,1.313) 


2u 3 4 5 

(0.1875,0.938,0.938) 


2u2 3 4 5 

(0.281,0.797,1.406) 


2v? 3u^ 4 5 

(0,2.064,0) 

18-4.2 

12 3 6 

(0.1875,0.75,1.875) 


lu 2 3 6 

(0.375,0.891,1.313) 


lu^ 2 3 6 

(0.281,1.172,1.031) 


lu 2u^ 3 6 

(0.5625,0.75,1.125) 


1 2u^ 3 6 

(0,1.875,0) 


1 2v? 3u 6 

(0.281,0.844,1.406) 


lu 2v? 3u 6 

(0.469,0.985,0.844) 


lu^ 2u^ 3u 6 

(0.656,0.422,1.406) 


lu 2u 3u 6 

(0.1875,1.5,0.75) 


1 2u 3u 6 

(0.1875,0.9375,1.313) 

18-4.3 

12 3 4 

(0.5625,0.9375,1.688) 


lu 2 3 4 

(0.281,1.781,0.844) 


lu 2u 3 4 

(0,2.625,0) 

18-4.4 

12 5 6 

(0.5625,0.9375,1.688) 


lu 2 5 6 

(0.281,1.781,0.844) 


lu 2u^ 5 6 

(0.75,1.125,0.75) 


lu 2u 5u 6 

(0.1875,1.6875,1.3125) 
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